arXiv: 1009.5327 



o 



Uniform Approximations for the M/G/l Queue with 
Subexponential Processing Times 

Mariana Olvera-Cravioto and Peter W. Glynn 

Abstract: This paper studies the asymptotic behavior of the steady-state waiting time, Woo, of the 
M/G/l queue with subcxponencntial processing times for different combinations of traffic intensities 
and overflow levels. In particular, we provide insights into the regions of large deviations where the 
so-called heavy traffic approximation and heavy tail asymptotic hold. For queues whose service time 
distribution decays slower than e _v ^ we identify a third region of asymptotics where neither the 
$-H heavy traffic nor the heavy tailed approximations are valid. These results are obtained by deriving 

approximations for P(W / C x> > x) that are either uniform in the traffic intensity as the tail value goes to 
infinity or uniform on the positive axis as the traffic intensity converges to one. Our approach makes 
clear the connection between the asymptotic behavior of the steady-state waiting time distribution 
' and that of an associated random walk. 
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1. Introduction 



We study in this paper the asymptotic behavior of the steady-state waiting time distribution of an M/G/l 
queue with subexponential service time distribution and first-in-first-out (FIFO) discipline. The goal is to 
provide expressions that will allow us to identify the different types of asymptotic behavior that the queue 
. experiences depending on different combinations of traffic intensity and overflow levels. We give our results 

' for the special case of an M/G/l queue with the idea that the insights that we obtain are applicable to more 

£C) . general queues and even to networks of queues. 

The special case of an M/G/l queue with regularly varying processing times was previously analyzed in [17], 
0^ ' where it was shown that the behavior of P(Woo > x), the steady-state waiting time distribution, can be fully 



described by the so-called heavy traffic approximation and heavy tail asymptotic (see Theorems 2.1 and 2.2 
in [17]). As pointed out in that work, the same type of results can be derived for a larger subclass of the 
subexponential family, in particular, for service time distributions whose tails decay slower than e _v/ *. As 
the main results of this paper show, the behavior of for lighter subexponential service time distributions 
may include a third region where neither the heavy traffic approximation nor the heavy tail asympotic arc 
valid, and where the higher order moments of the service time distribution start playing a role. The exact 
way in which these higher order moments appear in the distribution of is closely related to the large 
deviations behavior of an associated random walk and its corresponding Cramer series. 

The approach that we take to understand the asymptotics of P(W / o > x) over the entire line is to provide 
approximations that hold uniformly across all values of the traffic intensity for large values of the tail, or 
alternatively, uniformly across all tail values for traffic intensities close to one. From such uniform approxi- 
mations it is possible to compute the exact thresholds separating the different regions of deviations of , 
which for service time distributions decaying slower than e _v/ * are simply the heavy traffic and heavy tail 
regions, and, for lighter subexponential distributions, include a third region where neither the heavy traffic 
approximation nor the heavy tail asymptotic hold. Similar uniform approximations have been derived in 
the literature for the tail distribution of a random walk with subexponential increments in [6], [7], and [20], 
where the uniformity is on the number of summands for large values of the tail or across all tail values as 
the number of summands grows to infinity. The results in the paper are in some sense the equivalent for the 
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single-server queue. 

To explain the idea behind our main results let us recall that one can approximate the tail distribution of 
the steady-state waiting time of a single-server queue with subexponential processing times, P(Woo > x), 
via two well known approximations: the heavy traffic approximation and the heavy tail asymptotic 



respectively, where V\ denotes the service time, t\ the inter-arrival time, and p the traffic intensity of the 
queue. We refer the reader to Chapter X of [1] and the references therein for more details on the history 
and the exact formulation of these limit theorems. The heavy traffic approximation is valid for the general 
GI/GI/1 queue and can be derived by using a functional Central Limit Theorem type of analysis (see, e.g. 
[11, 12]). The theorem that justifies this approximation is obtained by taking the limit as the traffic intensity 
approaches one and is applicable for bounded values of x. The heavy tail asymptotic is valid for the GI/GI/1 
FIFO queue with subexponential service time distribution (see, e.g., [10]), and is obtained by taking the 
limit as x goes to infinity for a fixed traffic intensity, that is, it is applicable for large values of x. One can 
then think of combining these two approximations to obtain an expression that is uniformly valid on the 
entire positive axis. 

The approach we take in the derivation of the main theorems is to start with the Pollaczck-Khintchinc formula 
for the distribution of the steady-state waiting time of the M/G/l queue, which expresses it as a geometric 
random sum, and use the asymptotics for the tail distribution of the random walk. One of the difficulties 
in obtaining uniform asymptotics for the distribution of Woo hes in the highly complex asymptotic behavior 
of the random walk. Surprisingly, most of the cumbersome details of the asymptotics for the random walk 
disappear in the queue, but showing that this is indeed the case requires a considerable amount of work. The 
qualitative difference between queues with service time distributions with tails decaying slower than e _v/ * 
and their lighter-tailed counterparts comes from the asymptotic behavior of the random walk associated to 
the geometric random sum. The function e _v/ * has been identified as a threshold in the behavior of heavy 
tailed sums and queues in [6, 16], and [3, 13, 14], respectively, to name a few references, and we provide here 
yet another example. 

As mentioned before, the approximations we provide can be used to derive the exact regions where the heavy 
traffic and heavy tail approximations hold, but we do not provide the details in this paper since our focus is 
on deriving uniform expressions for P(W 00 > x) under minimal conditions on the service time distribution. 
The setting we consider is the same from [3, 14] where the busy period was analyzed. More detailed comments 
about the third region of asymptotics that arises when the service time distribution is lighter than e _v/ * can 
be found in Remark 2 right after Theorem 3.4. For clarity, we state all our assumptions and notation in the 
following section, and our main results in Section 3. 

Finally, we mention that the expressions given in the main theorems can be of practical use as numerical 
approximations for P{W ao > x), and based on simulation experiments done for service times with a Pareto 
(a > 3) or Weibull (0 < a < 1/2) distribution, they seem to perform very well (see Section 4 in [17]). It 
is worth pointing out that the uniform approximations given here are far superior than the heavy traffic or 
heavy tail approximations individually even in the regions where these are valid, which is to be expected 
since they are based on the entire Pollaczek-Khintchinc formula; they are also easy to compute given the 
integrated tail distribution of the processing times and its first few moments (cumulants). 

2. Model Description 

Let (W n (p) : n > 0) be the waiting time sequence for an M/G/l FIFO queue that is fed by a Poisson arrival 
process having arrival rate A = p/EVi and independent iid processing times (V n : n > 0). Provided that the 
traffic intensity p is smaller than one, we denote by Woo(p) the steady-state waiting time of the queue. We 
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assume that G(x) = P(V± < x) is such that its integrated tail distribution, given by F(x) = J* G{t)dt/ EV\ 
is subexponential, where G(t) = 1 — G(t). The sequence {^Q}i>i will denote iid random variables having 
distribution F. 

Define Q(t) = — log F(t) to be the cumulative hazard function of F and let q(t) = {EV\)~ 1 G(t)/ F{t) be its 
hazard rate function; note that q is the density of Q. Just as in [3] and [4], we define the hazard rate index 

r = limsup^. (2.1) 

All the results presented in this paper hold for subexponential distributions G (its corresponding integrated 
tail distribution F) satisfying the following assumption. 

Assumption 1. a.) < r < 1; 

'2, ifr = 0, 

4/(1 -r), if r ^0. 



b.) liminf t _ i . 00 tq(t) > a(r), where a(r) 



Assumption 1 is consistent with Conditions B and C in [3] and [4], respectively, and also very closely 
related to Definition 1 in [14] . All three of these works study the asymptotic behavior of random sums with 
subexponential increments applied to either the study of the busy period of a GI/GI/1 queue or to ruin 
probabilities in insurance. Also, by Proposition 3.7 in [3], Assumption 1 (a.) is equivalent to the function 
Q(t)/t r+s being decreasing on t > to > 1 f° r anv < S < 1 — r, which is the same as equation (3) in 
[20] , where uniform asymptotics for the tail behavior of a random walk with subexponential increments were 
derived. As mentioned in [3] and [4], Lemma 3.6 in [3] implies that sup{fc : Spf-f] < 00} > liminf^oo tq(t), 
so Assumption 1 (b.) guarantees that -E[Xf] < 00 for all k < a(r). Furthermore, Assumption 1 (b.) and 
Lemma 3.6 in [3] together imply that liminf t _ i . 00 Q(t) / \ogt > liminf t _>. 00 tq(t) > a(r), which in turn implies 
that for some (3 > a{r) > 2 and to > 1, 

Q{t)>p\ogt foralU>i - (2.2) 

Although the tail distribution of the busy period in queues with heavy tailed service times is related to that 
of its waiting time in the sense that it is determined by G{x) (see [3, 4, 14, 23]), the approach to its analysis 
is rather different from that of the waiting time, so the only connection between the results in this paper 
and those cited above is the setting. 

This family of distributions includes in particular all regularly varying distributions, F(x) = x~ a L(x) with 
a > 2, and all semiexponential distributions, F(x) = e~ x with < a < 1; in these definitions L 
is a slowly varying function. The regularly varying case with a > 1 was covered in detail in [17]. Some 
subexponential distributions that do not satisfy Assumption 1 are those decaying "almost" exponentially 
fast, e.g. F(x) = e~ x / l °s x . 

Before stating our main results in the following section, we introduce some more notation that will be used 
throughout the paper. Let fx = EX 1 = EV?/{2EV X ) and a 2 = Var(Xi) = EV?/{3EVi) - (EV? / (2EVx)) 2 . 
Also, define 

K = maxll e {0,1,...} : limsup ^ffi, > ol + 2, (2.3) 

t i->oo t l IV +L > J 

and note that by Proposition 3.7 in [3], Q{t)/t r+s is eventually decreasing for all 8 > 0, which implies that 
Q(t)/t r+s for all 6 > 0. In particular, for r £ [0, 1/2) this implies that Q(i)/t 1 / 2 -> and k = 2. Also, 
we obtain the relation (k — 2)/(k — 1) < r, or equivalently, n < (2 — r)/(l — r). Combining this observation 
with our previous remark about Assumption 1 (b.) gives that for < r < 1 and any 2 < s < (2 + r)/(l — r) 
we have E[X? +S ] < 00. 
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3. Main results 



As mentioned in the introduction, the idea of this paper is to use the Pollaczek-Khintchine formula to write 
the distribution of the steady-state waiting time as 

oo 

P(Woo(p) > x) = £(1 - P )p n P{S n > x), (3.1) 

71=0 

where S n = X\ + • • ■ + X n and {Xi}i>i is a sequence of iid random variables having distribution F, and 
then approximate P(S n > x) by an appropriate asymptotic expression. The approximation that we use for 
P{S n > x) depends on the pair (x, n), and for the heavy-tailed setting that we consider here, one can identify 
four different regions of deviations. 

The first region is the one described by the Central Limit Theorem (CLT), i.e., where 

P(S n >x)^l-^((x-nn)/y/ria), (3.2) 

and <£(•) is the standard normal distribution function. The second region is the so-called Cramer region, 
which provides additional correction terms to the CLT approximation. When the distribution F has finite 
exponential moments, the Cramer approximation is given by 

p((s n - > x) = exp , (E))( 1 + ( x '^+ l 



1 - $(a;/- v /n) V Vn 

where X(t) = J"^_ 3 XjP / j\ is a power series with coefficients depending on the cumulants of X\ known in 
the literature as the Cramer series (see [18], Chapter VIII, §2, or [15]). When F is heavy-tailed, nevertheless, 
X(t) diverges for all t and a truncated form of this series replaces A(-). In the setting of this paper, only the 
terms up to k (as defined by (2.3)) are needed, and we obtain the following approximation for P(S n > x) 

P(S n > x) ps (1 - $((£ - nfily/rur)) e K^) +"<2«( £ ^) i (3.3) 

where 

Q«(') = E^r. ( 3 - 4 ) 

A2 = —1, and {Xj}j>7i are the coefficients of the Cramer series corresponding to Y\ = {X\ — p)/a. Note that 
if k = 2, then approximations (3.2) and (3.3) are the same. 

The third region is known in the literature as the "intermediate domain" , and the exact asymptotics for 
P(S n > x) in this region can be considerably complicated (see [7] and [20] for more details). Fortunately, the 
range of values corresponding to this region in the Pollaczek-Khintchine formula is negligible with respect 
to the rest, and we will only need to use an upper bound for P(S n > x). The fourth and last region is the 
heavy-tailed region, also know as the "big jump domain" (sec [6] and [9], for example), where 

P(S n > x) ps nF(x — nfx). 

In the discussion above we purposefully omitted describing the boundaries between the four different regions, 
since that alone requires introducing various (complicated) functions and their corresponding asymptotic 
behavior. In terms of the Pollaczek-Khintchine formula, it is enough to consider simpler versions of those 
thresholds. We start by defining the functions 

ui(t) =t 2 /(Q(t)Vl) and u> 2 {t) = t 2 /{Q{t) V l) 2 , 

where a; V y — m&x{x,y} (x A y — min{a;, y}), and let w^ 1 (t) = inf{u > : t < u)i(u)}, i = 1,2. We give 
below some properties of the u>~ 1 operator; the proof is omitted but can be derived through straightforward 
analysis. 
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Lemma 3.1. For any continuous function uj : [0,oo) — > [0, oo) such that lim t _ i . 00 tu(t) = oo, define the 
function ui^ 1 : [0,oo) — > [0,oo) as = inf{w > : t < lu(u)}. Then, the following are true 

a. ) lu^ 1 is monotone non decreasing and left- continuous. 

b. ) uj^ 1 is a right inverse of ui , that is, a;(w~ 1 (t)) = t, for all t > w(0). 

c. ) iful{t) = sup 0<s<t uj(s), then uj^ 1 is a right inverse of ID for all t > lo(0). 

d. ) < t'for all t > 0. 



We now define the threshold functions delimiting the different regions of asymptotics for P(S n > x). Let 

K r (x) = 



l{x~u~\x))/p\V0, if re [0,1/2), 
\mm{u 2 {x) 1 x/(2p)}\ V 0, if r G [1/2, 1), 



M(x) = \{x~uj^ l {x))/p\ V0, and N(x) = [(x - ^x\ogx)/p\ V 0. 

Note that if r G [0, 1/2) and if 8 > is such that r+S < 1/2, then w 2 (t) > Ct 2(1 " r "' 5) for some constant C > 0, 
so Lo-\t) < Ct^W-r-W = o(t). Also, provided r + 5 e (0, 1), Wi(t) > Ct 2 ~ r ' s so ^(t) < Ci 1 '^^ = 
o(t). Therefore, all three functions are strictly positive for large enough x. Moreover, as mentioned in the 
previous section, Assumption 1 (b.) implies that Q(t) > f3\ogt for all t > t for some /3 > a(r) > 2, which 
gives u>2(t) < 0J\(t) < fi~ x t 2 1 logt, and ui^ix) > (x) > *J ((3/2)x\ogx for all x > Xo- We then have that 
K r (x) < M{x) < N(x) for all large enough x. 

To better understand the definitions of the threshold functions consider the zero mean case with finite 
variance, for which it is well known that the CLT approximation (3.2) holds for x = 0(y/n); translating into 
the positive mean case, this gives rise to the threshold n > (x~y/cx)/p for some constant c > 0. Substituting 
the constant by log x gives the threshold N{x). The Cramer approximation (3.3) holds, in the zero mean case, 
uniformly for x < o-\(n), where <7i(n) is the solution to the equation x 2 = nh(x) and E[e h ^ Xl ^ l(Xi > 0)] < oo 
(see, [8] §5.1 and the references therein); taking h = Q gives the threshold n > oji(x), and translating into 
the positive mean case gives n > (x — 1 (x/p))/p. Note that E[e^^ Xl ^ 1(X\ > 0)] = oo but, for example, 
E]fi Q{x)-2\ og Q(x) l( ^ Xi > o)] < oo, so this choice of h is very close to the boundary of the region. Finally, 
the asymptotic P(S n > x) ~ nF(x) as x — > oo is known to hold, in the mean zero case, for n < cw2(x) (see 
Theorem 1 in [4]), and provided that w^" 1 (x) = o(x) (which occurs when r e [0, 1/2)), the translation into the 
positive mean case gives the threshold n < (x — ui^ (x / p)) / p,. When r G [1/2, 1) we cannot guarantee that 
t02(x) < x/fx, so by taking the minimum between ui2(x) and x/(2p) we satisfy the condition n < loi(x — np), 
and therefore our choice of K r {x). We point out that since the thresholds do not need to be too precise, we 
ignored the constant p inside of wf 1 and uj^ 1 in the definitions of M(x) and K r (x), respectively, to simplify 
the expressions. 

The first asymptotic for P(W OD (p) > x) we propose is given by the following expression based on the 
Pollaczck-Khintchine formula, for k = 2, 

K r (x) 

Z K (p,x) = (! - p)p n nF{x -np) + E \p a ^ l(aZ < ^pw^(x)/^)] , (3.5) 



n=l 



and for n > 2, 



Kr(x) ,- N(x) nQ J 



Z K (p, a) = £ (1 - P)R n nF{x np) + ]T (1 



P)p n - 



2ttu ^— ' x — nu 

^ n=M(x) + l ' 



+ E \p a (*' z ) l( a Z < xZ/xloga;)] , (3.6) 



Y^ n =A a " — ® whenever B < A. Our first theorem is formally stated below. 



where Z ~ N(0,1) and a(x,z) = ix — az^Jxj ' p \ /p. Throughout the paper we use the convention that 
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Theorem 3.2. Suppose Assumption 1 is satisfied, and define Z K (p,x) according to (3.5) and (3.6). Then, 

P(W oc (p)>x) 



lim sup 

x-S-oo <p<l 



0. 



Remark 1 . (i) We point out that the approximation given by Z K (p, x) is explicit in the sense that given the 
exact form of F, all the functions and parameters involved in the approximation are known. In particular, 



E 



a(x,Z) 



\{oZ<JUT) 



= p» e 



^? $ [ xjz — | — Y— log p 



3/2 



(ii) This approximation is suitable for numerical computations since it involves no integrals or infinite sums. 

(iii) With some additional work once can show that the first term in (3.5) and (3.6) can be replaced by 

K r {x) 

F(x) £ (1 - p)p n n, 

n=l 

which is asymptotically equivalent to the heavy tail asymptotic pF(x)/(l — p) for appropriate values of (x,p). 
We choose not to use this simpler expression because our numerical experiments show that it would result in 
a less accurate approximation for P(W 00 (p) > x). (iv) For the case k > 2, the middle term in (3.6) provides 
a direct connection between the Cramer region of asymptotics for the random walk and the asymptotic 
behavior of the queue, and also reiterates the qualitative difference between distributions decaying slower 
than (re = 2) and those with lighter tails (see [16], [13], [14], to name some references), (v) Unlike the 

next approximation, given in Theorem 3.4, the expression Z K {p,x) does not work as a uniform asymptotic 
in x > as p /• 1 for P(Woc(p) > %), since it does not converge to one for small values of x. Nevertheless, 
it is not difficult to show that 

P{Woo(p) >X) 



lim sup 

P/* 1 x>x(p) 



for any x{p) — > oo as p 



Z K (p,x) 

1 (see the proof of Lemma 3.3 in 



1 







17]). 



In the same spirit of the heavy traffic approximations in [22] and [5], where P(W 00 (p) > x) is approximated by 
e xS(p) whgj-g S(p) is a power series in (1 — p), our second result derives an approximation that involves a power 
series in log p. The number of terms in this power series is also determined by n (as in the definition of Q K (-)), 
and its coefficients are closely related to those of the Cramer series. This other approximation substitutes the 
second term in (3.5) and the second and third terms in (3.6) by their corresponding asymptotic expression 
as p /* 1. The intuition behind this substitution is that these terms only dominate the behavior of Z K (p,x) 
when the effects of the heavy traffic are more important than those of the heavy tails. Besides unifying the 
cases k — 2 and k > 2, this new approximation will also have the advantage of being uniformly good for 
x > as p /* 1. In order to state our next theorem we need the following definitions. 

Let 

K p {t) = (1 - 1) log p + V V ^ ( ) ~ t\ (3.7) 



where A2 = —1, and {^j}j>3 are the coefficients of the Cramer series corresponding to Y = (X\ — p)/a. This 
function can be obtained by expanding (1 — t)Q K (po~~ 1 t/(l — t)) into powers of t; the details can be found 
in Lemma 6.1. We also need to define u(p) to be the smallest positive solution to A' p (t) = 0. Some properties 
of A p and u(p) are given in the following lemma. 

Lemma 3.3. Define A p according to (3.7) and let u{p) be the smallest positive solution to A' p (t) = 0. Then 
A p is concave in a neighborhood of the origin, 



n=l 
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as p / 1, where b\ 



A p (wO)) 
■ and for n > 2, 



logp+^(logp) 2 + 0(|logp| 3 ), k>2, 



dt n ~ 



t ~° (mi,...,m»-i)eA,-i 



/ 2 \ n + s n-i n— 1 -. 
£ (n + fln _i - !)!(-!)" ^ I — yfo l(j < K - 2))' 



Ai = {(mi, . . . , m n ) e N" : lmi + 2to2 + • • • + nm n = n} , s n = mi + • • • + m n , and 



K—2 



P K (t) = A' p (t) + log p±tJ2ajt j . 
The second approximation for P(Woo(p) > x) that we propose is 



A K (p,x) 



K r (x) 

Ed 



p)p n nF(x — np) + i 



(3.8) 



71=1 



where w(p, x) = min {u(p),uj^ 1 (x)/x}. The precise statement of our result is given below. 
Theorem 3.4. Suppose Assumption 1 is satisfied, and define A K (p,x) according to (3.8). Then, 

P(W 00 (p)>x) 



lim sup 

x ^°°0<p<l 



A K (p, x) 



0. 



Moreover, 



lim sup 



P(W 00 (p)>x) 



1 



= 0. 



A K (p,x) 

Remark 2. (i) As mentioned earlier, the difference between Z K (p,x) and A K (p,x) is in the terms that 
correspond to the behavior of the queue when the effects of the heavy traffic dominate those of the heavy 
tails. In particular, what prevents Z K (p,x) from being uniformly good for all values of x as p /*■ 1 is that if 
x is bounded, then the second term in (3.5) and the second and third terms in (3.6) do not converge to one 
when p /* 1, which can be fixed by substituting them by their asymptotic expression asp / 1; evaluating 
A p at the value w(p,x) = min{u(p),a; 1 " 1 (a;)/x} guarantees that the contribution of e» lAp ^ tu ' p '' r ^ becomes 
negligible when the queue is in the heavy tail regime, (ii) For analytical applications, Lemma 3.3 states that 
A p (u(p)) can be written as a power series in log/5 whose terms of order greater than k can be ignored. For 
numerical implementations, nonetheless, it might be easier to compute u(p) by directly optimizing A p (t), 
since A p (t) is just a polynomial of order k. (iii) By simply matching the leading exponents of the heavy tail 
asymptotic and the function -A p (u(p)), that is, by solving the equation 



/' 



logp = -Q(x) 



we obtain that the heavy tail region is roughly IZi = {(x,p) : p < e - ^^^" 1 }, whereas on IZ2 = {(x,p) : 
p > e~^^ x ^ x } one should use e^ Ap< -- u ^ to approximate P(Wao(p) > x). It follows that the heavy traffic 
region is given by the subset of IZ2 where e~ Ap< - tl< - p - ) ' ) is asymptotically equivalent to e~~^ 1 ~ p \ the heavy 
traffic approximation for the M/G/l queue. We note that when k = 2, the heavy traffic region is the entire 
IZ2, but it is a strict subset of 72.2 if k > 2, in which case a third region of asymptotics arises where neither 
the heavy traffic nor the heavy tail approximations are valid, (iv) As mentioned before, the coefficients of 
A p (t) can be easily obtained from the first k — 2 coefficients of the Cramer series of Y = (Xi — p)/a, which 
in turn can be obtained from the cumulants of Y. 



We end this section with a formula that can be used to compute the coefficients of the Cramer series. 
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3.1. Cramer Coefficients 

The following formula taken from [21] can be used to recursively compute the coefficients in the Cramer 
series, and we include it only for completeness. 

Proposition 3.5. Let Y be a random variable having EY = 0, Var(y) = 1, and cumulants 71,72,.... 
Let A3,A4, ... be the coefficients of the (formal) Cramer series of Y , i.e., X(t) = ^2jL 3 \jt J Let Aj = 
{(ni, . . . , nj) £ W : lni + 2n2 + • • • + jrij = j}. Then, for j > 3 and Sj-2 = fi\ + • • • + n,j-2, 



The rest of the paper consists mostly of the proofs of all the results in Section 3 and is organized as follows. 
Section 4 states an approximation for P(S n > x) that is valid for all pairs (x, n) and that will be used to 
derive uniform asymptotics for P(W oa (p) > x). Section 5 contains the proof of Theorem 3.2; and Section 6 
contains the proofs of Lemma 3.3 and Theorem 3.4. Wc conclude the paper by giving a couple of numerical 
examples comparing the two suggested approximations for the tail distribution of Woc(p), Z K (p,x) and 
A K (p, x), in Section 7. A table of notation is included at the end of the paper. 



4. Uniform asymptotics for P(S n > x) 

In this section we will state the uniform approximation for P(S n > x) that we will substitute in the Pollaczck- 
Khintchinc formula (3.1) outside of the heavy-tail region. This approximation was derived in [20] for mean 
zero and unit variance random walks and it works on the whole positive line as n — > 00. Although rather 
complicated as an approximation for P(S„ > x), it will be useful in the derivation of simpler expressions for 
the queue with the level of generality that we described in Section 2. For the heavy-tail region (small values 
of n) we will use in section 4.1 a result from [4] to prove that P(S n > x) = nF(x — np)(l + o(l)) as x — > 00 
uniformly in the region 1 < n < K r (x). 

We start by stating the assumptions needed for the mean zero and unit variance random walk, and after giving 
the approximation in this setting we will show that under Assumption 1, the random variable Y\ = (Xi—fx) /a 
satisfies these conditions. Then we will apply a slightly modified version of the approximation to the positive 
mean case and we will show that it holds uniformly in the region n > K r (x). 

The notation f(t) x g(t) as t 00 means < uminft_>oo f(t)/g(t) < limsup^^ f(t)/g{t) < 00. We will 

also use C to denote a generic positive constant, i.e., C — 2C, C = C + 1, etc. 

Assumption 2. Let Y be a random variable with E[Y] = 0, Var(y) = 1 and tail distribution 




The first four coefficients arc given by 



A3 = 73, A 4 = 7 4 - 37 3 , A 5 
A6 = 76 - 157573 - IO74 + 1057473 



75 - IO7473 + 157!, 
105 73 4 



1 - V(t) = V(t) 



t —> OO, 



where D(t) = J , <t t 2 



dV(dt), Q has Lebesgue density q, and satisfies 




t->oo Q(t) 



and 



liminf Q{t)/\ogt > f/(l - f). 
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Suppose further that £ , [|F| K+1 ] < oo, where 

K = max I I G {0, 1, 2, . . .} : limsup ^U+D > f + 2 - 

Throughout this section let Qa(t) = YTj=2 \^ 7i'> wnere ^2 = — 1 and {Xj}j>3 are the coefheients of the 
Cramer series of Y, and let S n — Y\ + ■ ■ ■ + Y n , where {F,} are iid with common distribution V(t). We also 
define the functions 

b{t) = t 2 /(Q(t) V 1), and b'^t) = in£{u > : t < b(u)}. (4.1) 
We start by proving some properties about the functions Q, and o _1 . 

Lemma 4.1. Suppose Assumption 2 holds. Then, for any s € (f, 1) there exists a constant to > 1 sitc/i i/iai 

a J Q(t)/t s is decreasing for all t >to, 

b. ) b-\t) < iV(2-«) / or allt>t , 

c. ) & _1 (ci) < c 1/(2_s) 6 _1 (t) < c6 _1 (t) /or aZZ i > i a™d any c > 1, 

d. j 6 _1 (ci) > c6 _1 (t) /or t >to and any c < 1, 

Afeo, i/ie following limit holds 

e. ) Um t ^ 00 e-'5( t / fe " 1 W)g(6- 1 (t)) = 0, 

Proof. Part (a.) follows directly from Proposition 3.7 in [3]. For part (b.) note that Q{t)/t s is eventually 
decreasing for any f < s' < s, so 

hm < sup — hm ; = 0. 

It follows that Q(t) < t s for all t > to for some to > 0. This in turn implies that b(t) > t 2 ~ s for all t > to, 
and therefore, b' 1 ^) < iVC 2 "*). 

For part (c.) note that Proposition 3.7 in [3] gives Q(co _1 (i)) < c s Q(o _1 (i)) for any c > 1 and all sufficiently 
large t, then 

V V " Q(cV(2- S )6-i(t)) " Q(6- x (*)) 

It follows from noting that b(t) is strictly increasing for large enough t, that 

c l/(2- 8 ) 6 -l (t) > 

For part (d.) let c < 1 and define u(x) = c _1 o(x), u(x) = o(c _1 x). By Proposition 3.7 in [3], Q(x) > 
c"Q(c x), from where we obtain that 

, , c~ 1 x 2 c(c~ 1 x) 2 (c~ 1 x) 2 , . 

u(x) = -= < — ^ '- = i '- < v(x). 

Q(x)Vl ~ c s Q(c- 1 x)Vl c- 1 + s Q(c~ 1 x)V c- 1 ~ 

It follows that > v~ 1 (x), where u~ 1 (x) = mi{t > : cx < b(t)} = b~ 1 (cx) and v~ 1 (x) = 

inf {t > : x < b(c~H)} = cinf{t > : x < b{t)} = cb~ l {x). 
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For part (e.) let v = liminf t _ i . 00 Q(t)/\ogt > f/(l — f) and note that 

k-l/_A\2 



lim e-^/^^Qib-^t)) = lim e-QW^ffl) (b W) = Urn e-OW^Wi/r'W) g (*)) 



U 2 



lim e -QW«)/«) " = i im e -«(«/«(«»Q( u ) 

0(u) u-foo 



< lim e- l/log(u/( 3 (tl)) Q(M) 

u— >-oo 

lim 



Q(«) 



jf/H-i) I 

By part (a) Q(u) < Cu s for any s > f and w sufficiently large, and by assumption f < vj(y + f), so simply 
choose f < s < vj{y + 1) to see that the last limit is zero. □ 



Lemma 4.2. Suppose Assumption 2 holds. Define 

L(h)= e ht dV(t), and H(z) = inf (nln£(/i) - zh). 
J-oo h>0 

Then, for any constant c > 0, 

e H(z) =e nQ k (f) (1 + (1)) (4.2) 

as n — > oo, uniformly for y/n < z < c£> _1 (n). 

Proof. Choose < <5 < i — f and set s = f + 5. Define r](z) = b~ 1 (z 2 ) and 

n(z,n) = (i - $(z/y/n)) l{z < y/n) + (l - *(z/y/n)) e^ +H ^ . (4.3) 

Suppose first that f G [0, 1/2) and note that in this case k = 2 and nQ^(z/n) = — z 2 /(2n). Note that we can 
choose 6 above so that s < 1/2. Then, by Lemma 4.1 (a.), Q{t)/t s decreases for all sufficiently large t. Also, 

V(—z) = for all z > fi, 



D(z) 
and 

/n/n(y/n) roc 
u 2 dV{t) =1-1 _ («(«) + 2/M)e- Q ( tl) dw 
-A* Jn/ri(y/n) 

= 1 + (e-'S^MVS)^ =l + (l/Q(r7(V^))) (by Lemma 4.1 (e.)) 
as n — > oo. Define Xn = & _1 (n) = ij{y/n) and note that 

Xr 2 t = HXn) = 1 

Q(Xn)n n 

Then, by Lemma la in [20], we have 

n(z,n) = {l-*{z/y/n)) (1 + o(l)) 

= (1 - $(z/y/n)) e^+" Qs (^)(l + o(l)) 
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as n — > oo, uniformly for ^Jn < z < "fXn, where 7 > is an arbitrary constant (see the statement of Remark 1 
in [19] to see that the constant 7 can be arbitrary). 

Suppose now that f <G [1/2, 1) and recall that by assumption £?[|y| K+1 ] < 00. Then, by Lemma lb in [20], 

n(z,n) = (l-$(z/V^))e E "=i^W^ (l + o(l)) 

= (1 - $(z/y/n)) e^+" Qs (™)(l + o(l)) 

asn-> 00, uniformly for y/n < z < jr](y/n), where 7 > is an arbitrary constant. To see that 7 can be 
arbitrary see Remark 1 in [19] where the statement of the result is 

e ff(*) =e nQ„(*) (1 + (1)) 

as n — > 00, uniformly for < z < A„, for a function A„ that in [20] is taken to be A n = r](y/n), and verify 
that all the arguments go through if we let A„ = r){y/yn) for any constant 7 > 0. Then, use Lemma 4.1 (c.) 
to see that 77(^/7") = b~ l {"jn) < (7 V l) 1 ^ 2 -^-*^). □ 

The main approximation is given below. 

Theorem 4.3. Suppose Assumption 2 holds. Fix e £ (0, 1) and set 

ir(y, n)=(l- HvIM) 1(V < Vn) + (1 - *(v/Vn)) e^ +nQ ^) l(y > Vn), (4.4) 

J(y,n) = Vn~l v(t)^' 1^1) dt+ -= V{t)e nQ ^) dt \, (4.5) 



3/-V" V V ?1 / V 27T i v Sv( r l)- 1 (2(l+e)ii)) 

C„ = mm t \\ + ^n) . (4.6) 

t>v^r y2 i 2 J 

Then, as n — > oo, uniformly in y, 

P (S n > y) = (Hy, n) l(y < (1 + e)C„) + J(y, n) l(y > (1 - e)C n )) (1 + o(l)). 
Moreover, there exist constants < 71 < 1 < 72 such that C n £ [71& (n), 7 2 6 _1 (n)]. 

Proof. Choose < 6 < 1 — r and set s = r + 5. Note that by Lemma 4.1 (a.) Q(t)/t s is eventually decreasing. 
Also, since Var(Yi) = 1, 

P (§ n < ty/inj -> *(*) 

by the CLT. Define L(h) and H(z) as in Lemma 4.2 and let ir(z,n) be given by (4.3). 

Set r](z) = b^ 1 (z 2 ) and note that by Lemma 4.1 (b.) b~ x (t) < t 1 ^ 2 ^^ for all t sufficiently large, so r/(z) = 
o(z 2 ). Since D(t) — > 1 as t — > 00, we have 

D{z 2 /r](zj) = 1 + o(l) = D(z), z^oo. 

Let 7 = 1/(2(1 + e)), and define 

u n = b n/7), v n = sVnQiVn) . 

L>(V^) 

Then, by Theorem 2 and Remark 1 from [20], 

P (s n > y) = (*(y, n) i(y < (1 + e)C„) + { jf ' V{t)& dt 

1 rV—Vn ~) \ 

7(4)^* 1(2/ > (1 - e)C n ) (1 + o(l)) 



27T J^\f{ y -\) 
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asn-> oo, uniformly for all y and for any A € v n \. Also, by Lemma 4.2, 

e ff(») = e nQ *(#)(l + o(l)) 

uniformly for s/ri < z < ci> _1 (n) for any c > 0. We will show below that C„ < & _1 (2n) < 2 1 /( 2 ~ s '6~ 1 (n) 
(by Lemma 4.1 (c.)), so we can replace ?r(y, n) by 7f(y, n). Also, by choosing A = 6 (2(1 + e)n) and noting 
that for t > y — A we have y — t < A = 6 _1 (2(1 + e)n) < 4 1 /( 2_s )& _1 (n) (by Lemma 4.1 (c.)), we can replace 
e H(y-t) yfftfa e "Qs( jL s-)_ This gives the statement of the theorem. 

To verify the order of magnitude of C n let h(t) = t + ^r-nj and note that h is continuous and a.s. 

diffcrcntiablc. Recall that by assumption Q has Lebesgue density q, and note that b(t) is eventually increasing, 
since by Lemma 4.1 (a.) Q(t)/t s is eventually decreasing. Then, for all to < t < 6 (2n(l — s)), 

1 Q(t) - tq(t) 1 1 

h (t) = n ■ s < nil — s) ■ —rr < 0. 

w 2 i 2 - 2 v ; b(i) - 

For t > b~ 1 (2n) note that liminf t _ i . 00 tq(t) > liminf t _ i . 00 atq(at + fi) — 2 > a(r) — 2 > 0. It follows that 

1 Q(t) - 1 1 

fi (t) = — — n ■ - > — — n • — 7-T- > 0. 

w 2 i 2 ~ 2 6(i) ~ 

We conclude that C n € 1 (2(1 — s)n), & _1 (2n)], and by by Lemma 4.1 (c.) and (d.), 

r 1 (2(l-s)n)>(2(l-s)Al)r 1 (ti) and 6" 1 (2n) < (2 V l) 1/(2_a) b _1 (n). 

□ 

We now give a lemma stating that under Assumption 1, the random variable Y± = (Xi — \x)jo satisfies 
Assumption 2. Throughout the rest of the paper, 

Q{t) = Q{at + fi)-2logt, (4.7) 

and the functions b and b~ l , as well as the constant f, are defined according to this function. 
Lemma 4.4. Suppose Q satisfies Assumption 1, then Y\ = (X\ — fi)ja satisfies Assumption 2. 

Proof. Let Q(t) = Q(at+fi)-2logt, thenF(i) = P(Yi > t) = er^/t 2 , and since D(t) = f H<t t 2 dV{t) -> 1 

as t — > 00, then V(t) X D{t)t~ 2 e~^ t \ Also, since Q has Lebesgue density q, then Q has Lebesgue density 
q(t) = aq(at + fi) — 2/t. It follows that 

- ,. tq(t) taq(<rt + fi) - 2 zq(z) 
r = lim sup — — = urn sup — < lim sup 



t^oo Q{t) t-yoo Q(crt + fi) - 2 log t z^oo Q(z)-21ogz 

< r Qjz) 

< r lim sup — — — . 

z^oo Q(z) -2 log z 

By (2.2), there exists f3 > a(r) > 2 such that Q(t) > j3 log t for all sufficiently large t. It follows that 

r lim sup — < 



1 -2(log z)/Q(z) 1-2//3' 
where if r > we have r/(l - 2/fJ) < r/(l - 2/a{r)) = 2r/(l + r) < 1. Therefore, r < f < 1 and 

liminf^ = liminf^±^-2>^-2. 

t-»oo logt i->oo logt 
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Clearly, if r = then f = and /3 — 2 > = f/(l — f). If r = we already showed that f < /3r/ (/3 — 2), which 
combined with fj > a(r) = 4/(1 — r) gives f < 1 — 2/(/3 — 2), which in turn implies that /? — 2 > 2/(1 — r) > 
f/(l - r). 

We also note that for any i G {0, 1, 2, . . . } 

lim sup ® (t) - a 1 '^ lim sup ~ 2 log " + 2 logg 

Since for any / € {1, 2, 3, . . . } we have limsup^^ Q{t)/t l ^ l+1 ^ = a l ^ l+1 ^ limsup.^^ Q(u)/u l ^ l+1 \ it fol- 
lows that R = k. Finally, from the discussion following the definition of k, equation (2.3), we have that 
E[X? +S ] < oo for any 2 < s < (2 + r)/(l - r), which implies £[|Yi|" +1 ] < oo. □ 



We are now ready to give a uniform approximation for P(S n > x) that will work over the region n > K r (x). 
We choose not to use this approximation in the heavy tail region 1 < n < K r (x) to avoid having to show 
that it is equivalent to the heavy tail asymptotic nF(x — nfi). Instead, we use a result from [3] that will give 
us without much additional work the heavy tail asymptotic directly. 

We point out that we will not apply Theorem 4.3 to the positive mean exactly the way it is stated, but 
instead we use a slight modification that will work better when applied to the queue. In particular, we will 
substitute the function 7r(j/,n) given by (4.4), where y = (x — njS)/a, with the following 

f$(—y/y/x/(i), k = 2, 

7t K (x, n) = | H _ y/ ^ } 1(B > N{x)) + ^ e ^(f) l (n < N(x)), k > 2. ^ 



The function J(y,n) given in (4.5) does not need to be modified since its contribution will be shown to be 
negligible in the queue. 

Lemma 4.5. Suppose Q satisfies Assumption 1. Let y = [x — n[i)/a, fix e € (0, 1) and define 

B K (x, n) = tt k (x, n) l(y < (1 + e)C„) + J(y, n) l(y > (1 - e)C„), 
where 7r K (x,n), J(y,n) and C n are given by (4.8), (4.5) and (4.6), respectively. Then, 

P(S„ > x) 



lim sup 

x ^°°n>K r (x) 



0. 



B K (x,n) 

Moreover, there exist constants < 71 < 1 < 72 such that C n £ [71b -1 (/in), 726 _1 (/xn)]. 

Proof. By Theorem 4.3 and Lemma 4.4, we have that 

P(S n >x) = (7r(y, n) l(y < (1 + e)C n ) + J(y, n) l(y > (1 - e)C„)) (1 + o(l)) 

as 1 -> 00 for all n > K r (x), where Tf(y,n) is given in (4.4). Furthermore, by the same theorem and 
Lemma 4.1 (c.) and (d.), there exist constants < 71 < 1 < 72 such that C n <S [jib^ 1 (/in), ^b -1 (nn)]. It 
can be verified that 

{y < (1 + e)C n } C{y< 2 72 6- 1 (/in)} C {x - 2a l2 b- 1 (x) < n/i} = {n > l(x)} 

for sufficiently large x, where l{x) = (x — 2aj2b~ 1 (x)) / ft, so all that remains to show is that fr(y,n) = 
7r K (x, n)(l + o(l)) asi->oo for all n > l(x). 

Note that after some algebra we can obtain the equivalence 

X 0~ ' \J~X 

{y < Vn} = {n > m(x)}, where m(x) = - + 
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Since N(x) = [(x — yfx log x) / fi\ < m(x) for sufficiently large x, it follows that for k > 2, 
\% K (x,n) - 7r(y,n)| l(n > l(x)) 

= <E> (^-y/x/x/fij l(n > N(x)) - $ (—y/y/n) l(n > m(a;)) 

g ™Q*(*) < n < N(x)) - $(-y/V^)e^ +nQ ' t (") l(Z(x) < n < m(x)) 



< 



$ {—y/y/xfjA - $ {—y/y/n) l(n > m(x)) 
$ (-y/y/x/fij - $(-y/y/n)e£ +nQk (%) l(JV(a;) < n < m(a;)) 

e »^(f ) _ ^{-y/ y /T 1 )e^ +nQ ^) l(l(x) <n< N(x)), 



yy/WfJ, 

while for k = 2 we have y 2 /(2n) + nQ K (y/n) = and 

|* re (x,n) -7f(y,n)|l(n > = $ ( -y/y/xjj?) - $> (-y/y/n) l(n>l(x)). 

To analyze (4.9) and the corresponding segment of (4.12) define s(x) = (x + -y/xlogcc)/^, then 



(-y/y/x/jl) - $ (-y/y/n) l(n>m(i)) 

$ (y/y/x/fj) - $ (y/y/n) l(m(x) < n < s(x)) + 2$ (y/y/n) l(n > s(x)) 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



<*'(o)M 



l(m(x) < n < s(x)) + 2$ f— (s(x)fi - x)/y/a 2 s(x)j l(n > s(x)) 



C (a: „T )2 l(m(z) < n < s(.t)) + 2$ f ) l(n > «(a:)) 



,3/2 



a(l + o(l)) 



< Cmin { $ (- ^ SX ) } 1(„ > 

Since for n > m(x) we have $ (—y/y/x/Ji) > $ f— (a; — fim(x))/y/ a 2 x/ ftj — > $(—1), it follows that (4.9) 
and the corresponding segment of (4.12) are bounded by 

CVi(a;)$ (^-y/y/x/Jlj l(n > m(x)), 

where </?i(:c) = min {(log x) 2 /y/x, $ (— v / /71oga;/(2cr)) }. To bound (4.10) and the corresponding segment 
of (4.12) we note that for N(x) < n < m(x) we have nQ K (y/n) = —y 2 /(2n) + 0(y 3 /n 2 ) (recall that 
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-Xy/n) = -2/7(2") if/s = 2),so 



1-e' 



V / Jn.-r.ln. n z 



y/nx/fj, 



y/y/x/n 



< C 



where for the third inequality we used the relation $(— z) > ^'(z)z/(z 2 + 1) for all z > 0. Therefore, (4.10) 
and the corresponding segment of (4.12) are bounded by 

Cip 2 (x)$ (—y/y/x/Jn l(N(x) <n< m(x)), 

where <P2(x) = Qogx) 3 / 2 /y/x. To bound the last segment of (4.12) note that the preceding calculation yields 

$ \-y/\/x/Jx\ - $ (-y/y/n) l{l(x) <n< N(x)) 



< C 



< c 



(y 2 + x) (x — n/j.) 



$ -y/y/x/fi l(l(x) <n< N(x)) 



{{b-^x)) 2 +x)b- 1 (x) 



$ {-y/\fx/Jj\ l(l(x) <n< N(x)). 



Since n = 2 implies that Q(t)/y/i — > 0, then 



lim 



(b-\x)f 



x—¥oo X 



lim 

X— ¥ OO 



{b-\x)f 



= lim 



so 



where 



C 



b{b~ x {x)) t^oc t 2 /Q(t) 
((b-Hx)) 2 +x)b^(x) ^Ab~\x)Y 



lim f3(x) = lim 



(b-\x)f 



< c- 



lim 



lim^=0. 



(6(6-i(a;)))2 t'-T™ (t 2 /Q(t)) 2 t^Zo t 
We have thus shown that when k = 2, 

$ (-y/y/x/n) - $ (-y/\/n) l(n > i(s)) < C max ^(a;)$ (-y/y/x/n) l(n > Z(x)) 
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Finally, to bound (4.11) we use the inequalitvy §'(z)z/(z 2 + 1) < <&(— z) to obtain, for l(x) < n < N(x), 



< 



< 



1 



< C 

< c 



V\fU y 2 /n+lj V2tt 

n 3/2 

V x/u — y/n H 7T- 

x — nfi x 



yV27r 



x (x — nfi) 2 J y^fhFfl 



1 



log x _/ y^fhijl 



It follows that (4.11) is bounded by 



C Vi (x)^= e nQ ^) l(l(x) <n< N(x)), 



where ipi(x) = b~ 1 (x)/x + 1/logx. We conclude that 



\tt k (x, n) — 7r(x, n)\ < C max tpAx)Tr K (x, n). 

ie{l,2,3,4} 



for all n > l(x). This completes the proof. 



□ 



4-1- A first approximation for P(W oc ,(p) > x) 

We will now give an approximation for P{W ao {p) > x), that although too complicated to be used in practice, 
will serve as an intermediate step towards obtaining the more explicit approximations given in Theorems 3.2 
and 3.4. 

The idea of this section is to substitute P(S n > x) in the Pollaczek-Khintchine formula (3.1) the heavy-tail 
approximation nF(x — nfx) in the range 1 < n < K r (x), and by B K (x, n), as defined in Lemma 4.5, in the 
range n > K r (x). 

The intermediate approximation for P(W 00 (p) > x) is given by 

K r (x) oo 

S K (p,x) = J2 (l-p)p n nF(x-nri+ £ (1 - p)p n n k (x,n)l(y < (1 + e)C n ) 

n=l n=K r (x)+l 
oo 

+ (l-p)p n J(y,n)l(y>(l-e)C n ), (4.13) 

n=K r (x) + l 

where y = (x — pn)/a, and n K (x, n), J{y,n) and C n are given by (4.8), (4.5) and (4.6), respectively. The 
last term in (4.13) corresponds to the so-called "intermediate domain", where as mentioned in Section 3, the 
asymptotic behavior of P(S n > x) is rather complicated. Under additional (differentiability) assumptions on 
Q, more explicit asymptotics for J(y,n) have been derived in [20] (see also [7] for other results applicable 
to this region). We point out that S K (p,x) is "very close" to being the approximation in Theorem 3.2 if we 
replace l(y < (1 + e)C n ) with l(n > M(x)) and ignore the entire third term of S K (p,x), to see this sum the 
tail of the second term of S K (p, x) to write it as the expectation of a function of a normal random variable. 

We will now show the asymptotic equivalence of P(Woo(p) > x) and S K (p,x). 
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Lemma 4.6. Suppose Q satisfies Assumption 1, then 

P{S n > x) 



lira sup 

l<n<K r (x) 



nF(x — nfj,) 



- 1 



= 0. 



Proof. Recall that u>2(x) = x 2 /(Q(x) V l) 2 and w 2 (x) = inf{w > : x < ^(u)}- By Lemma 3.1, u> 2 1 is non 
decreasing, W2(w^~ (x)) = x and w^" 1 (o;2(a;)) < x. Let t n = w 2 _1 ((/iA l)n)/2, and note that 



hm sup v n — - — i < lim sup v n — . ^— - — - < 



2 /-QCwa 1 ^)) 
— —77- lim sup v s = 

w 2 x ((mA l)n)/2 " (fiAl) 1/2 s ^oo ^ \s) 



(MAI) 



-—lim sup 



{w 2 (w 2 " 1 ( S ))} 1 /2 ( M A 1)V2 



Then by Theorem 3.1 in [4], 



lim sup 



P{S n -fxn>t) 



nF{t) 



- 1 



= 0. 



Next, we will show that for n < K r (x) we have x — fin > t n . 

First, when < r < 1/2 we have K r {x) = [(x — ui^ 1 (x)) / fi\ , so n < K r (x) implies 

X — fM > X — flK r (x) > U>2 (x) > W 2 ~ 1 (/ X ^r(^)) > tR r (x) > ^n- 

Similarly, when 1/2 < r < 1 and K r (x) = [min{w2(^), x /(2fi)}\ , we have that n < i^ r (a;) implies 

x — fin > x — fiK r (x) > max{x — jjLU2{x), x/2} > a; 2 ^ 1 (a;2(a;))/2 
> Lu 2 ; 1 (K r (x))/2 > t KAx) > t n . 



These observations, combined with the fact that the subexponentiality of F implies that P(S n > x) = 
nF(x)(l + o(l)) as x — > 00 uniformly for 1 < n < a(x) for some a(x) — > 00 completes the proof. □ 

Combining Lemmas 4.6 and 4.5 gives the following result. 

Proposition 4.7. Define S K (p,x) according to (4.13) and suppose Q satisfies Assumption 1, then, 

P(W 00 (p)>x) 1 =Q 
S K (p, x) 



lim sup 

1 - >00 0<p<l 



This first approximation for P{W ao (p) > x) might not very useful in practice since it involves two integrals, 
those in the definition of J(y, n), that are not in general closed- form, and two indicator functions that depend 
on the quantity C n (the solution to a certain optimization problem). The approximation given in Theorem 3.2 
is more explicit, and thus more suitable for computations, both numerical and analytical. 



5. Proof of Theorem 3.2 



The proof of Theorem 3.2 is rather technical, so we divide into several lemmas, the first of which gives some 
more properties of the functions 6 _1 and w^f 1 . 

Lemma 5.1. Suppose Q satisfies Assumption 1. Let Q and b~ x be defined according to (4.7) and (4.1), 
respectively. Then, 
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a. ) lim^oo Q(i)/wr 1 (t) = lim^ Q{t)/b-\t) = 0, 

b. ) lim^oo Vt/cu^it) = 0. 

Proof. To show the first limit in (a.) use Proposition 3.7 in [3] with some r < s < 1 as follows, 

Q(t) 



lim 

t—^-oo 



< ( U /9(«))'QM = lim fQ(«) 



it— >-oo \ i/, 



0. 



For the second limit we first note that the same arguments used above give limt_>. 00 Q(t)/b *(0 = 0, so all 
we need to show is that limsup^^ Q(i)/Q(t) < oo. That this is the case follows from 



< hm sup — < lim sup - , 

" t^oo Q(t) ~ *->«, Q(ot)- 2 log t 



Q(t) ... (a' 1 V lYQju) C 
< hm sup = lim sup 



Q(u) - 21ogu + 2 logcr 

and (2.2), which gives Q(u)/ log u > (3 > a(r) > 2 for large u. 
For part (b.) 



l-2bg«/Q(«) 



lim ; = lim = = lim = lim 



(*) 



(*) 



ti— >-oo ^ 



o. 



□ 



Next define Z K (p,x) according to (3.5) and (3.6), and S K (p,x) according to (4.13). Let 
Ei(p,x) -- 



Ei{p,x) 
Ez{p,x) 



J2 (l-p)p n K(x,n){l(y<(l + e)C n )-l(n>M(x))} 

n=K r (x) + l 

E„°Lm W+ i(1 - pK^M - E [p<^ 1 (aZ < ^\x)/^)] 
ZZn^+A 1 - P)P n *«M - E [p<*' z ) 1 (aZ < VTHc^)] 

OO 

£ (l-p)p n J(y,n)l(y>(l-e)C n ). 

n=K r (x)+l 



re = 2, 
re > 2, 



Then, 



1 5« (p, x) - Z K (p, x) | < Ei (p, x) + E 2 (p, x) + E 3 (p, x) . 

We will split the proof of Theorem 3.2 into three propositions, each of them showing that Ei(p,x) = 
o(Z K (p, x)) as x — y oo uniformly for < p < 1, and some auxiliary lemmas. We start by giving a result that 
provides lower bounds for Z K (p, x). 



Lemma 5.2. Fix c > and let p(x) = e C ^Q( X )/ X . Then, for any < p < p(x), 



C p 

Z K (p,x) > F(x), 

1 - P 



w/w7e for p(x) < p < 1, 



^(p^O^Ce-^W")), 
where A p (u(p)) was defined in Lemma 3.3. 
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Proof. Let J{x) = \ xj \J Q{x)\ < K r {x) and note that 

J{x) 

Z K (p,x) > ^(1- p)p n nF{x-np) 

n=l 

J(x) 

> (1 - p)F(x) J2 np n 



= F(x)-^—{l-p J ^-(l-p)J(x)p 
The first statement follows from the observation that for < p < p(x) we have 



p m < 



>-J(x) _ e -c^Q(x)+o(l) _^ g 



For the second statement consider first the case K = 2, for which e^VWPJJ = e f"*n—& an d 
Z K (p,x) > E \p a ^^ l(aZ < yfiiu?(x)/y/x) 



e ■ log p+ $ f v^r^) + ^| log p 



crWx 



p: 



3/2 



> e tA P («(p))$ f vm^eI _ = e tM«(P)) (1 + o(1)) 

\ o^fx y/p% I 



as x —> oo, for all < p < 1 (since Q(x)/lo 1 1 (x) ~ > by Lemma 5.1 (a.)). For k > 2 we split the interval 
1) into two parts as follows. Define p(x) = e a ' 2 ^ . Then, for p(x) V p(x) < p < 1, 

Z k ( P ,.t) > £ [p<^> z ) 1 (Z < VJIbi^/tr)] > $(0)e* log ' ,+ ^ r:B(log '' )2 > CetM«(P)). 
For the interval [p(x), p(x) Vp(x)) (assuming p(x) > p(x)), let it n = (x — np)/x and use Lemma 6.1 to obtain 



/-vi ^ -a„(m„) 

G{l z_ p) j2 >cq 



g w (p,g) > ~ v >: — — >c(i-p)Vi 

' u N(x)-l 

By Lemma 3.3, A p is concave on [0,um(i)], and its maximizer, u(p), satisfies 



«M(x) „fA p (u) 



«(p) = -- logp + 0(1 logp| 2 ), A p («(p)) = logp + ^ (logp) 2 + 0(| logp| 3 ). 
Also, the derivatives of A p satisfy 

2 2 

A'„(i) = -logp- ^t + 0(t 2 ) and A"(t) = ~ + 0(t). 
Then, for some £t between i and u(p) and some constant C, > p/a 2 , 

A p (t) = A p («(p)) + ^^(t - U (p)) 2 > A„(u(p)) - ^(t - U (p)) 2 . 
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Note that for p < p(x) we have u(p) > \f\og x/x + O (log x/x) . Therefore, for any < 8 < 1 and x sufficiently 
large, 

(1 — p)\fx I du 

> (1 - p)Vie* A '< u <"» / 6 rfz 

> Ce tA P («(p)) ^ _ $ (-^^(p) - U^j.!))) 

> Cei Ap{u( - p)) (since Va?u(p) ^ oo for p < p(x)) . 

□ 

The next lemma will be useful in showing the uniformity in < p < 1 of our bounds. 

Lemma 5.3. Let a(n,x) be any function that does not depend on p. Then, for any l(x) > 4p~ 1 x/Q(x) and 
m{x) < ^, we have 

m(x) m(x) 

/or sufficiently large x. 

Proof. Define p(x) = e~^^ x ^ x . By Lemma 5.2, we have that for < p < p{x), 

i(x) m(x) 



'g (l-p)p»a(n,x)< sup C 'g'a-^V" 1 ^'^ 



o< P <p(x)^(p,a;) n f^ x) o<p< P (,) n=l{x) F(x) ' 

Define /i„( i o) = (1 - p) 2 p n ~ l and compute h' n {p) = (1 - p)p n ~ 2 (n(l — p) - I - p). Note that for p e (0, /5(a;)] 
we have 

n(l- p)-l- p> n(l - /3(x)) - 1 - p(x), 

so > on (0,p(x)] for all n > (1 + p(x))/(l - p{x)) (note that (1+ p(x))/(l - p{x)) ~ 2p- 1 x/Q(x) as 

x — > oo). Therefore, 

m(:c) . . m(x) 

sup £ (l- p )V- 1 %^= £ (l-p(x)) 2 / 5( a ;r- 1 e^)a(n,x) 

< CQ [ X) Y e^-^a(n,x). 

n—l{x) 



M. Olvera-Cravioto and P. Glynn /Uniform Approximations for the M/G/l Queue 



21 



For the range p(x) < p < 1 fix e S (0, 1) and use Lemma 5.2 again to obtain 

m(x) m(x) 

sup t~7 T H (l-p)p n a(n,a:) < sup — g2 V {1 - p)p n a{n,x) 

p(.)<p<l 2„(P, X) p{x) < p<1 e * logp+Cl-^^MloSP) 2 

< y, SU P e M SF a(n,a;) 

= o(^xj ^ ^(SH^^a^a:) (for all n < x/n) 

n—l(x) 
r<n>t \ m ( x ) 

CQ(x) ^ (£_„)£2M , s 
= > ' x am, a;). 

x ^ 

n—l(x) 

□ 

Proposition 5.4. Under the assumptions of Theorem 3.2, 

r Ei(p,x) 
lim sup 75—7 r- = 0. 

x^oc Q<p<1 Z K {p,X) 

Proof. Define m e (x) = min{n G {1,2,...} : np+(l + e)aC n > a;}, and recall that M(x) — [(x — oj^ 1 (x))//j,\ . 
Let y = (x — np)/a. Then, 

00 

Ei(p, x) = (1 — p)p n it K (x, n) l(min{TO c (x), M(x)} < n < ma,x{m e (x) , M (x)}) . 

n=K r (x) + l 

Choose < 5 < 1. By Lemma 4.5 there exist constants < 71 < 1 < 72 such that C„ e \lib~ 1 (pn), j2b~ 1 (pn)] . 
Then, for any n > l(x) = (x — r )\ab~ 1 (x)) / p and x sufficiently large, 



V < lib-^x) <(-Vl) -fib~\pn) < (1 + e)C n , 
\pn ) 



where in the second inequality we used Lemma 4.1 (a). Similarly, for any n 
and x sufficiently large, 

y > 2l2b^{x) > 2 l2 b- 1 {pn) > (1 + e)C„. 
It follows that L^( a; )J < m c (x) < [l(x)\ for sufficiently large x. Hence, 

max) [l(x)\,M(x)} 

Ex(p,x)< (1 - p)p n Tt K (x,n), 

n=min{ [k(x)\ ,AI(x)} + l 

and by Lemma 5.3, 

I \ \ mnx{[l(x)],M(x)} 

Ei{p,x) CQ(x) ^ Q , x) _ ynQM . , , 

sup — -< > e WK ' ir k {x,n). 

0<p<l A K (p,x) x . f-f 

n— nim{ [k(x)] ,M{x)\ + l 

By using the inequality <fr(—z) < $'(z)/z for any z > 0, and observing that n = (x/p)(l + o(l)) for all 
min{ [fc(^)J , M(x)} <n< max{ [l(x)\ , M(x)}, we obtain, for such n and all sufficiently large x, 

< _>^L e -# + -^Le"^(£) < C^-e-^ 1 ^ < Ce^ 1 " 5 ^. 
y^l-np yy/Zirp y^fp 
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It follows that 



SU p E W' X ) < °Q( X > y e ^M-(i-«)^ 

o< P <i Z K (p,x) ~ x . . fr 1 . ... x 

n— mm{ \k\x)\ ,M(a:)} + l 



. . max{ [/(a:)J ,Af (a:)} 2 / 2<t 2 Q(k) 

CQ(X) s—y ~( 1 ~ l5 ) i 2V l^ 1 " |l-Jl,i(.^(l(.)vM(.)]) 



X J 

ra=min{Lfe(x)J,M(aj)}+l 



Note that by Lemma 5.1 (a. 



l im ®W < c lim 9M_ = 0, (5.1) 

x^oox- n([l(x)}V M(x)) ~ x^oo Ao^x) 



which implies that for sufficiently large x, 



d / ^ l~ir\l \ ma,x{[l(x)],M(x)} 

Ei(p,x) CQ(x) ^ 

SU P V - e 

0<()<1 ^ ft 1 z . f-^ 

n=min{ L«;(x)J ,M + l 

CQ(x) f»-Mmin{Lfc(x)J,M(x)} _ 2J£ ^ 

< / e ^ ; so** <xu 



^ ii-fimax) LZ(a;)J ,M(a;)} 

< ^ f e"^. 

J (1 ;^/ ?r ( 7lt Tb-i(^)A^- 1 (^)) 

Finally, by using the inequality $(— z) < $'(z)/z for z > again, and (5.1), we obtain that 

Um sup — < lim — i = 0. 

x^oc 0<p<1 Z K (p,x) x^oo 'f 1 ab^ 1 (x) Auj 1 (x) 



□ 



Lemma 5.5. Let h K (x) = 1 (x) if K = 2 anii h K (x) = \fx log x, if K > 2, i/ien 



>!C )<^V^p(«-M S 



/i K (x 



Proof. Recall that a(x, z) = // 1 fx — azy/xj /ij and Z ~ N(0,1). Define L K (x) — \_(x — h K (x))/p\. Note 
that exact computation gives, 



£ (l - p)p n $ {-y/V^ 

=L K (x) + l 



= E 

= E 

= E 

= E 



J2 (1 - P)p n HZ > y/^/xTp) 

n=L n {x)+l 

max{ \a{x,Z)\+l,L K ,{x)+l} 



p La(x,z)j+i ! ^ Z < ( x _ ^(x))/^^)] + (-(x - pL K (x))/^/^xJpj 
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Observe that h K (x)/ 'y ' a 2 xj 1 p < (x — \iL K (x) ) / y/a 2 x/p < h K (x) / yj a 2 x/ p + p?l 2 / %/ cr 2 x, from where it follows 
that E2 (p, x) can further be bounded by 



E 2 (p, x) < E f p La(»,2)J+l 1 (h K (x)/y/^xJJl < Z < (x - (iL K (x)) / y/^xjjl 



E [( P wa)j+i _ 1 ( z < ^/V^Z") 



Next, note that since a(x, z) is decreasing in z, we obtain that (5.2) is bounded by 



.3/2 



\fx P 



For (5.3) we use the simple bound 



(1 - p)£ [p a ^ 1 (Z < h K {x)/^x~/J? } 
And for (5.4) we use the inequality $(— 2;) < $'(z)/z for any z > to obtain the bound 



\ / h K (x) 



e!M«ir 



(5.2) 
(5.3) 
(5.4) 



□ 



Proposition 5.6. Under the assumptions of Theorem 3.2, 

E 2 (p,x) 
hm sup — = 0. 

x^ca Q<p<1 Z K {p, X) 



Proof. Let h K (x) = lo 1 (x) if k = 2 and h K (x) = sfx logx, if ft > 2, then, by Lemma 5.5, we have that 



h K (x) 



(1 - p) E p<*^ 1 (aZ < ^JIh K (x)/^) 



Fix c > 1 and define /5(x) = e 



We will first show that E%(p, x) is o(Zk(p, x)) as x — > 00 uniformly for 
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< p < [>{%)■ Before we proceed note that h K (x)/ \J~x — > oo as x — > oo (by Lemma 5.1 (b.) for k = 2) and 



E 



p< x ' Z H(aZ <^JJh K (x)/Vx) 



< e" 



\ a^fx p 3 / 2 



+ e S lo - + ^^lf|logp|<^M + 41 

where for the inequality we used §{—z) < $>'(z)/z for z > 0. Furthermore, 

e" of+ ^ 1 log/3 < - W + ^— ■ p= < e A 1 2- 
\ cHa; fry a; / 

It follows that for sufficiently large x, E2(p,x) is bounded by 

Now we use Lemma 5.2 and the observation that h K (x)/x — > as x — > oo to obtain 

E 2 (p,x) p (x-h K (x))/n + J^-Mr 1 -^)/^ 

sup -^-7 r < C sup _ 

o<p<p(2:) o<p<pOr) P(1-PJ e w( - ' 

<Ce*W sup ^-^W-^/^+p^-^-^-^)/^ 

0<p<p(.T) V / 

= CeQW (e-^-^W-^^ +e -i( x - hs ^ 1 -^-^ BB ^ 21 
< c ( e -{^-^)Qi*) + e-i^-^-^^A _> o 



as x — > oo. 

For the range p(x) < p < 1 wc first note that Z K (p,x) > E [p a ( x > z )l UrZ < y/p,h K (x)/y/x)~\, so we have 

(1 - p)E [p a ^ 1 (aZ < ^ph K (x)/^)] 
sup L - —) — v - < 1 - p(x) -> 

p(x)<p<l ^k(Pi x ) 

as x — > oo. To analyze the remaining term we use Lemmas 5.2 and 3.3 to obtain Z K (p,x) > C~ Ap<,u< < p ^ > 
c , e *(u«/H-^ r (i«wrt a -ui(«>a)|iQg/,| 8 ) for gomc 7] > It follows that 

sup 



p(x)<p<i h K (x) Z K (p,x) 
< e 2<r^x sup 



|5(a)<P<l e f (logP+^i(logp) 2 -')l(«>2)|logp| 3 ) 

— i— e 2<^ x sup e" u ^ » . (5.5) 

h K \Z) 0<s<cpQ(x)/x 
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When k = 2 and h K (x) = lo 1 1 (x), (5.5) becomes 



-1/ ^\ 2 



sup e " va! / < 



Wf (x) 0<s<c M Q(a;)/x (x) 



which by Lemma 5.1 (b.), converges to zero as x — > oo. When k > 2 and h K (x) = \fx log x we split the 
suprcmum and bound (5.5) with 

C n log g I \/g Jog g s ct 2 x g 2 , 7?X „3 

e 2°' A < sup e fl f ' ' 



Vlog^ i^O<s<min{c/iQ(a')/a:^ _1/3 } 



sup e " 2^ f l(x /J < cfxQ(x)/x) 

x~ 1 / 3 <s<c[iQ(x) fx 

<^L supe"H^ + sup e-^+^^^l 1 ^ 

Vlogx 



s>0 s>a:- 1 /3 



2^ ^2^,1/6 



, _C I -f£(i 
< — == < 1 + sup e 



VlOgX [ 8 >x-V3 

as x — » oo. This completes the proof. □ 



Lemma 5.7. Let y = (x — nfj,)/a and fix < (5 < 1/2 and < c < 1. Define eg = S 1 (4/i 1 V l)er. Then, 
under the assumptions of Theorem 3.2, for all n < (x — cb" 1 (x)) / n and x sufficiently large, 

J(y,n) < Cn~F{ay + M )e^ Q(<ry) + ^—F(aVTi + ^e"^ 1 '^ l(un > x - c^iT^x)). 

2/ 

Proof. Let = &" 1 (2(1 + e)t/fi) and F(i) = P(Xj > at + fx) = F{at + n). Then from (4.5) we obtain 

that 

fl _ A(y-^/E)Af! e (un))/^/E _ , s 

J{y, n) < nV(y) +n V(y — sfrz)®' (z) dz + Cn \ V(y - zy/n)e nQh \ ^> dz. 



To analyze the integral involving Q K first note that if k = 2, then nQ K (z/ \/n) = —z 2 /2, while if k > 2 then 
nQ K {z/\ 
we have 



nQ K {z/\/n) = — z 2 /2 + 0(z 3 / \fn). Therefore, for 1 < z < ((y — y/n) A f3 e (iu,n)) / y/n and n sufficiently large 



nQ K {z/y/n) < -i 



from where it follows that 



'•((y-v / n)Aft(Ai™))/v / "_ (l-5)* 2 

J(y,n) < nV(y) + Cn I V(y - z\fn)e 2 dz. 



o 



We now bound the remaining integral with 



Cn I V(y- z\/n)e-^-^-dz (5.6) 



Cn / V(y-zVn)e ^ dz l(/3 e (/m) > 5y) (5.7) 

JSy/^/n 
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We start by analyzing (5.7), which is further bounded by 

CnV(y/n) / e 2 — dz l(/3 e (/m) > Sy) 

JSy/^n 



< CnV(V^)$ (- " /T Jf V ^j m{x) > Sy) 

, -,/n (l-S)S 2 y 2 

< CnV(V^) — e~ l(3 € (x) > Sy) 

V 



Cn 3 ' 2 ^ 



y 



V(y/n)e ~e 2 " V 1 (i-w.J l(/3 e ( x ) > Sy), 



where in the second inequality we used the relation z) < $>'(z)/z for z > 0. To obtain the second term 
in the statement of the lemma note that for n < (x — cfe~ 1 (x))/// we have 

2anQ(x) 2a 2 Q(x) 



(1 - S)S 2 xy ~ (1 - tf^Vcft-^a:) ' 
which converges to zero as x — > oo by Lemma 5.1 (a.). Then, for sufficiently large x, 

(l-6)6 2 y 2 I 2*nQ(g) \ (1 _ 5 )2 5 2 2 j2 2 

e 2„ \} (i-i)i^J < g- 1 — k < e —&- (S<l/2). 

Also, by Lemma 4.1 (c), /3 e (x) < (4/_i _1 V l)6 _1 (a;) = a~ 1 5c s b~ 1 (x). If follows that (5.7) is bounded by 

v(\m)e x Sn lum > x — cgb (x)). 

y 

To bound (5.6) we first note that by Assumption 1, q(t) < (r + S)Q(t)/t for sufficiently large t. Also, by 
Proposition 3.7 in [3], Q(t)/t is eventually decreasing, so we obtain 

V(y -u) = V(y)eW*% + . ^ dt < Ffoje^**^". 
Then, the change of variables u = z^/n yields the bound 

pSy 2 

C^V{y) I e SZ + -% + .^ e - il ^du 







< C^V{y) / e^^-^d, 
Jo 



/ - / , r\ Q(&y) (l-5)u- 

<CV^V(y) e {r+S) T±m u -^- du 



o 



= CW(y)e / e -~dz. (5.8) 

Now, define the set A = |cri/ > ui^ 1 ( ^rz^p^j- J } and note that t 2 /Q(t) = u±(t) is eventually increasing. It 
follows that for large enough x, 

f (r+j)V^ l f (r+jfQW . Sn 

A Q r (ffy) - (i - m~» / = 1 fw - ~ 

Also, 

, cr f , (r+^xl f (r+j)nQ(gj/) ^(1-*)^ 
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We then have that for z(x, n) = ^=2 { ^\^^y V) - Sy] , (5.8) is bounded by 
CnV(y) { e »■ 1(A) + e W^Fv* ~ sup l(A c ) } 

t>z(x,n) * (t) J 

< CW(y) { e—^r- 1(A) + e — ^ ^ — sup -^-f 1(A C ) 



Finally, we note that on A c = > x — 1 f ^r^p^H }, an d for sufficiently large x, (l-S)nfJ,/ ((r+S) 
1, so (5.8) is bounded by 

CnV(y) le— 1(A) + e" sup -^-^ 1(A C )| < Cny(y)e = . 



x) > 



□ 



Proposition 5.8. Under the assumptions of Theorem 3.2, 

E 3 (p,x) 
hm sup — = 0. 

x^oo Q<p<1 Z K (p,X) 

Proof. Set y = (x — nfx)/ 'a and recall that by assumption there exists ft > a(r) > 2 such that Q(t) > f3\ogt 
for all sufficiently large t. Now choose < S < min{(l— f)/2, (1 — r — 2//5)/2}. Note that by Lemma 4.5 there 
exists a constant < 71 < 1 such that C n > 7i&~ 1 (/in). Define c e = (1 — e) 2 7i<7 and c,5 = <5 _1 (4/i _1 V l)cr. 
Then, for any n > l e (x) = (x — c e b~ 1 (x))/ ft and x sufficiently large, 

V < (1 - e) 2 7i& -1 (z) < f — V lVl - e) 2 ^" 1 ^) < (1 - e )C„. 
Therefore, {y > (1 - e)C„} C {n < / e (x)}, and 

Li e (x)j 

n=K r {x) + l 

By Lemma 5.3, 

E 3 {p,x) CQ(x) L ^ )J £»oiEi 
SU P V- 2^ e * J (^> n )> ( 5 - 9 ) 

for sufficiently large x. Define ms(x) = (x — csb~ 1 (x))/ p. Then, by Lemma 5.7, 

Q£) ^ e ^l J(yin) <^) J- e ^M nF(CT2/ + /i)e ^ Q(CT!/) (51Q) 

ra=KV(x)+l n=K" r (a;)+l 



CQ(x) L '^ )J n 3 / 2 ^ ,V , 

2 F(CTVn + ^)e"~. (5.11) 



x e — ' y 

n=[m s (x)\+l 
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We start by showing that (5.10) converges to zero. To do so we first bound it with the following integral 

CQ(x) r KAX) \^F{u)e J ^ 1 ^du 

< CQ{x) / e ^-QM+^QMdu. (5.12) 

Now, by Proposition 3.7 in [3] we have that Q(x) < (x/u) r+s Q(u) for all u < x, from where it follows that 
uQ(x) nil ^ | 5(x — u) rt/ ^ ^ nil x ( (u\ x ~ r - & ^ S(x — u) 



~Q(u)+ S -^±Q(u)<Q(u)((™) 
x \\x/ 



<Q(u)((l-r-6)(^-^ ■ 5{X - U) 



= _Q(u)(x-u) {l _ r _ 25) 

X 

Let r\ = Next we will split (5.12) into three integrals and use one of the above inequalities to 

bound the exponent as follows 

r min{x-uK r (x), V x} /, _ , 5<m-u)\ 

CQ(x) / e QW((f) ~ 1+ ^)du (5.13) 

J X — fll e (x) 

+ CQ( X ) / e ~^^(l-r-28) du (514) 

J min{x — fiK r (x) ,rjx} 
px—fiK r (x) 

+ CQ(x) / e-^^^^-^du. (5.15) 

J min{x — fiK r (x) ,x/2 J rfi} 

To see that (5.13) converges to zero we note that it is bounded by 

p7]X POD 

CQ(x) e-Q^^du < CQ(x) e'^-^^du 

J X — fJ,le.(x) J X — fll e (x) 

= CQ(x)(x~[il c (x)r^ 1 - 2 V +1 

< CQ [ X } . (since p{\ - 26) > 2) 
x — fil e (x) 

CQ(x) 

where the last expression converges to zero by Lemma 5.1 (a.). To see that (5.14) converges to zero note 
that it is bounded by 



CQ(x) / e x^-^du 1(.t - fiK r (x) > rjx) 

J T}X 

CxQjx) r^d-W«W 
= — —. — r- / e dv l(x — fiK r (x) > rjx) 

Qyqx) J_( 1 _ r _25)(i_ T ,)Q( r , K ) 

CxQjx) liz^l QM Ll^m log(vx) 

~ Q(vx) ~ Cxe 

where in the last inequality we used Proposition 3.7 in [3] to obtain Q(x) < rj~ ( - r+s '>Q(rix) and then the 
assumption Q(t) > /3 log t. The last thing to notice is that our choice of 5 guarantees that (1 — r— 26)0/2 > 1. 
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Next, to analyze (5.15) we follow a similar approach and use the fact that Q{t)/t is eventually decreasing to 
obtain the bound 

AX—^K r (x) _ 

CQ{x) / e ? {1 - r - 2S) du l(x - pK r (x) > x/2 + fx) 

Jx/2+fj, 

< CQ(x) / e 3 {1 - r - 2S) du l{fx{K r {x) + 1) < x/2) 

J x/2 



Cx 2 



< ^^ e - " Q( 1f- w (i-r-2^) i( M ( Xr(a;) + i) < a./2). (5.16) 
i£ r (x) 

Now note that fx(K r (x) + 1) < x/2 implies that r € [1/2, 1), since for r <G (0, 1/2) we have x — pK r (x) = o(x); 
and in this case, 

{fi(K r (x) + 1) < x/2} C {jnm{(LL} 2 {x),x/2} < x/2} = {puj 2 (x) < x/2}. 
It follows that (5.16) is bounded by 

l ( 1 -'- 2 *) 1(huj 2 {x) < x/2) < CQ{x) 2 e-^ (1 ~ r " 2S) 



Cx yUMaaM /i „ 

g 2~ 

w 2 (x) 



< Cx 2 e 5 * o, 

where in the second inequality we used Proposition 3.7 in [3] to obtain that Q(x) < x r+s for large enough x. 
Finally, to prove that (5.11) converges to zero we first bound it with 

CQ(x) /-(^HDm ,3/2 s2 )2 

-t (cry sj fi)e ds 



*^ Jms(x)fi x s 



X-{l e (x) + l)n V / Jm s {x)n 



CxQ{x)—( , — -—\ 
F[a^{x)) 



6-!(x) 



oo 



< , , ,\ F (uJmsix)] I e ~ 2 dz. 



&JTt(x-(h(x) + l)u)/(2a^) 



Clearly, the last integral is bounded by a constant, and for the other terms we have 

lim xF (<jy/ms(x)) = lim t 2 F{t) = 0, 
since -E[X 2 ] < oo, and, by Lemma 5.1 (a.), lim. ! ;_ J . 00 Q(x)/b^ 1 (x) = 0. This completes the proof. □ 



Proof of Theorem 3.2. Propositions 5.4, 5.6 and 5.4 give 

S K (p, x) 



lim sup 

x-s-oo <p<l 



- 1 



= 0, 



which combined with Proposition 4.7 give 



lim sup 

x-s-oo 0<p<1 



P(W p (oo) > X 



Z K (p,x) 



- 1 



= 0. 



□ 
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6. Proof of Theorem 3.4 

In this section we prove Lemma 3.3 and Theorem 3.4. To ease the reading we restate the definition of A K (p, x) 
below. 

K T {X) 

A K (p, x)= (! - P)p n nF(x - np) + e i K " {wM) , 

n=l 

where A p is given by (3.7), w(p,x) = mhi^u(p),u>^ (x)/x} and u(p) is the smallest positive solution to 
A' p (t) = 0. 

We start with the proof of Lemma 3.3 and then split the proof of Theorem 3.4 into three parts. 
Proof of Lemma 3.3. That A p is concave in a neighborhood of the origin follows from 

2 2 

A' p (t) = -logp- ^t + 0{t 2 ) and A£(t) = + 0(t). 

If k = 2 we have A p (t) = (1 — t)\ogp — ^pit 2 , which is maximized at u(p) = —j^logp and satisfies 
A P («(p)) = logp+^(logp) 2 . 

In general, for n > 2 recall that P K (t) = A' (t) + logp, so u(p) is the solution to the equation P K {t) = log p. 
By Lagrange's inversion theorem, 

jn-i / + \™ (log/))™ 



M=0 



where 



i=2 j=2 J v J/ j=0 



Furthermore, by Faa di Bruno's formula, 
d 71 - 1 ( t 



b n = 



dt"- 1 \P K (t) / 

n — 1 1 



(mi ro n _i)£A„_i 



where _4j = {(m l7 . . . , ntj) G N-* : Lti! + 2m 2 + • • • + jrrij = j}, Sj = m\ + • • • + rrij, and ao = — ^j. Note 
that 6i = — a 2 /p 2 . Finally, since 

A p (i) = (l-i)log P -^ + 0(t 3 ), 



we have 



A p (u(p)) = log P + ^2 (log P) 2 + 0(1 log P\ 3 )- 



□ 



We now prove two preliminary results before we proceed to the proof of Theorem 3.4. 



M. Olvera-Cravioto and P. Glynn /Uniform Approximations for the M/G/l Queue 



31 



Lemma 6.1. Let A p be given by (3.7) and set u n = [x — nfj,)/x. Then, 



as x — > oo. 



sup 

0<p<l 



EN(x) n e 

n=M(x)+\ P x-ull 



N{x) 



- 1 



-M{x) + 1 xu n 



Proof. Define the function 



and note that 



A p (t) = (l-t) logp+(l-t)Q„ 



N(x) 



fit 



a(l-t) 



a; — n/u 



X ' ' U r , 

n=M(x) + l 



n=M(x)+l 

By expanding 1/(1 — t) J into its Taylor series centered at zero we obtain 



(l-t)Q K 



fit 



a(l-t) 



1 



j=2 ■ ?! °' J i=0 



j=2 i=0 

EE9 



x ^ ( r - i y+o{t^) 



j\ai \r - j 



EE, 



x ^ 3 ( r - i y+ (t^). 



r=2 j=2 



j\ai \r - j 



Recall from Section 2 (after equation (2.3)) that k < (2 — r)/(l — r) and Wj (a;) < Ca; 1 /' 2 r 5 ) for any 
< S < (1 — r) 2 /(2 — r) and a: sufficiently large. It follows that for < t < Um(x)+i —■ ^1 (%)/ x we have 



xt k+1 <xf 



: + l 



as x — } oo. Hence, 



Ar(z) 



E p n - 



n=M(x) + l 

as x — > oo, uniformly in < p < 1. 



< Ca; 



(a; 2 - 



•«) \ T=^T -(l-r) 2 +<i(2-r) 

= CX C2-r-i)(l-r) _> 



_ £3 Li P \ ' 



E 



n=M(x) + l 



-(1 + 0(1)) 



□ 



The second preliminary result is an application of Laplace's method, which states that the asymptotic 
behavior of an integral of the form 



-x<j>(t) 



f(t)dt, 



as x — y oo, is determined by the value of the integral in a small interval around the maximizer of 4> on the 
interval [c, d]. What makes the proof below very technical is that the limits of integration are functions of x. 
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Lemma 6.2. Let p(x) = e c > 0, k > 2, and define 7(x, p) = \J p log x/a + a^/xlog pj ' p?/ 2 . Then, 

under the assumptions of Theorem 3.4, as x — > oo, 



sup 



N(x) 

E 

n=A/(an) + l 



#1 s f 1 -*" 



x — np 



e fAp(«(p))$(_ 7 (p ;X )) 



-> 0. 



Proof. Let u„ = (x — np)/x and define c(x) = \/\og x/x, d{x) — oj 1 1 (x)/x. Then by Lemma 6.1, 



N(x) 



^ E p - 

n=M(x) + l 



iQ K (i-Z£ii) 



x — np 



1 iV(x) £A.(u„) 

1 ■r—v p ^ ' 



E 



n=M(x) + l 



(1+0(1)) 



' n=Af(x)+l ■ yu «+i 



■£ft(l + 0(l)) 



M y c (a;) * 

as i-> oo, uniformly for p(x) < p < 1. It remains to show that 



-dfc(l + o(l)) 



sup — 

p<p<i ei' 



A p (w(p)) 



(7v^(i - p) r w e 



rfi _ e t A p("(p))$(- 7 ( P;X )) 



c(x) 



We start by computing the derivatives of A p (i): 



A^) = -logp-4* + 0(t 2 ) and A^i) = + 0(t), 

and note that < -> for all < t < Also, we have u(p) = -^logp + 0((logp) 2 ) = - p) + 

0((1 — p) 2 ) e (0,d(x)) for all p(x) < p < 1. Set e = e(x) = 1/ log log x and note that for p(x) < p < 1, 
w(p) = o(d(x)), so for sufficiently large x we have 



£(P,*) = 



<j\fx(\ — p) /" txJ e 



p 3 / 2 \/2lT J c ( x ) t 



-d£-etM«W)$(- 7 (p,a:)) 



< 



oVx(l - p) 



c(x)v(l+e)«(p) e ^A p (t) 



/i 3 / 2 V27T ic(x)V(l-e)«(p) * 



dt- e ^ A " (tl(p)) $(- 7 (p,x)) 



Cv^(l-p) 



d(a) 



JA P (i) 



'c(a;)V(l+e)«(p) * 

To bound (6.1) note that for some £t between t and tt(p), 



c(^)V(l-e)u(p) e fA p (t) 



dt 



(6.1) 
(6.2) 



c(a:) 



c(x)v(l+e)u(p) e fA p (t) 



-dt — c 



A P (u(p)) 



c(a;)V(l-e)ii(p) 



o(x)V(l+e) U (p) g^fi^Ct-^Cp)) 2 



c(»v(l-e)u(p) 



so (6.1) is bounded by e^ Ap(u(p)) F(p, x), where 
F(p,x) 



3 / 2 V2^ 



c(z)V(l-e)w(p) 



dt - $(- 7 (p,x)) 
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To sec that suppr x \ <p<1 F(p,x) —> note that for (1 — s)u(p) < t < (1 + e)u(p) we have Ap(t) = — ^ + 
0(| logp|) and also t = tt(p)(l+o(l)) - logp(l+o(l)) = £(l-p)(l+o(l)). Let A = {c(x) < (l-e)u(p)) 
and let £(p) = max (i-e)ti(p)<t<(i+£)«(p) |A"(t)|. We start by analyzing F(p,x) 1(A), for which we have 

F(p,x)l(A) 

fxp, 



< 



(l + o(l)) 



e -fc-(t-«O0) dt _i 

(l-e)u(p) 



<<l- 



(i + o(i)) 



'"Hp) i \ ,2 



-eu(p) 



dz + $ 



+ $( 7 (p,x))|l(A) 

7tx(l — e)u(p) (Jsfxlogp 
a u, 3 / 2 



1(A) 



<{l- (l + o(l)) (l-2$ (-^W 
(1 + 0(1))") ">0 



_« (1 + OKp)/£)) K 1(A) 



< 3$ 



pxec(x) 



as x — > oo, uniformly for p(x) < p < 1. To analyze F(p, x) 1(A C ) we note that on A c we have e ^ ^ u( - p ^ 
e -^(*-(p)) 2 +o(- 2 «(p) 3 ) = (1 + (i))e-^( t -«W) 2 , which yields 

F(p,x)l(A c ) 
Ixp 



< 



(l + o(l)) 
/2x 



(l+o(l)) 



c(z)V(l+£)u(p) 



e 2o 



r(*-« (p)) (ft_$(- 7 (p,x)) 



c(x) 
F(cW-uW) 



1(A C ) 



e 2 



dz l(c(x) < (1 + e)u(p)) - §(-7(p, a)) 



l(A c 



<I> 



(u(p)-c(x)) -$(- 7 (p,x)) 



l(cW<(l + £)u(p))l(#) 



ipx 



+ $( _vj_ £u ( p ) ) i( c (a;) < (l + £ ) w (p)) + $(- 7 (p,x))l(c(x) > (l + e)uO»)) + o(l) 



< C 



(u(p) - c(x)) +7(p,x) 



l(c(x)<(l + e)u(p))l(A c ) 



* - 



£c(x) ) + $ I - 



/IX 



a(l 



■ec(x) + 0(xc(x) 2 ) + o(l) 



ec(x)(l + o(l)) + (1) ->0 



a(l + e 
< C^^x) 2 + 2$ 

V ° 

asi-} oo. We have thus shown that (6.1) i s o ^e~ Ap< '"' p ' ) ' ) ^ as x — > oo, uniformly for p(x) < p 
need to show that the same is true of (6.2). 
Note that (6.2) is bounded by 

r d(x) 



< 1. We now 



C7i(l - p) fet Ap(c(x)v(l+e)u(p)) + et A p ((l- £ )„(p)) < ^ _ / * I 

^ ' Jc(x) t 



(It 



<cyi(i-p)io g 



^1 Jgj ^ /' P tA p ((l+e) tl (p)) 

\/xlogx 



l((l + e)«(p)>c(x)) 



+e^ 



A P (c(x)) j^j + £)u( ^ < ^ + et A p ((l-e)«(p)) j^j _ £)u( ^ > c ^ 
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Note that for any t and some £t between t and u(p), 



A p (t) = A p (u(p)) - |A;(t)||«(p) - t\ + WM(u(p) - t) 2 , 
which gives that for (1 ± e)u(p) > c(x), A' p ((l ± e)u(p)) = T 2 ^^ + °( u (p) 2 )> and 

A p ((l ± e )«(p)) < A p (u(p)) ~ |A;((1 ± e)u{p))\eu(p) + ^e 2 u(p) 2 (l + (1)) 



= A p (u(p))-^s 2 u(p) 2 (l + o(l)) 
<A p (u(p))~^e 2 c(x) 2 (l + o(l)). 



For (1 + e)u(p) < c(x), A' p (c(x)) = -^(c(x) - u(p)) + 0(c(x) 2 ), and 

f> 2 



A p (c(x)) < A p (u(p)) - \A' p (c(x))\\u(p) c(x)\ + ^(u(p) c(x)) 2 (l + o(l)) 



A p (u(p))-^(c(x)-u(p)) 2 (l + o(l)) 
<A p (u(p))-^e 2 c(x) 2 (l + o(l)). 
Therefore, (6.2) is bounded by 

CV5(1 - ^(togsJe^GOJ-^c^i+oCi)) = o (et^P)) 

This completes the proof. 

Finally, we give below the proof of Theorem 3.4. 

Proof of Theorem 3.4-. Note that by Theorem 3.2 we know that 

0, 



lim sup 

x^oo o<p<l 



P(W p (oo) > x) i 
Z K (p,x) 



so for the first statement of Theorem 3.4 it suffices to show that 

A K (p,x 



lim sup 

a;^oo 0<p<1 



Z K (p, x) 



1 



0. 



□ 



(6.3) 



The second statement, which refers to the uniformity in x as p / 1 will follow from Lemma 3.3 in [17] once 
we show that supo <a!< (i_p)-i/4 \A K (p,x) — 1| — > as p /* 1. To see this is the case simply note that for all 

< x < (1-p)- 1 / 4 



K r (x) 



\A K (p,x)-l\<(l-p) J2 



\h p (w{p,x)) 



- 1 



< C(l - p)K r (x) 2 + 



< C(l - p)x 2 + Cx\ logp| < C(l - p) 1 ' 2 + (7(1 - p)- 1 ' A \\ogp\ ->■ 
as p 1. We now proceed to establish (6.3). 
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Let h R (x) = ui 1 (x) if k = 2 and h K (x) = y/x log x if n > 2, and set u n = (x — np)/x. Then, 



\Z K (p,x) - A K (p,x)\ 



N(x) 



y/'2irfi 



E ( 1 ~p)p- 



n=M{x) + l 



X — Tip 



1(k > 2) 



A p (u)(p,x)) 



< 



oyS(i - p) 

y/2irp 



N{x) 



E ^ 



n=M(x)+l 



g } e fA P («„) 

n=M(x)+l XU ™ 



1(k > 2) 



n ^ -A„(«„) 

e ^^i(k>2) 

v ^ n=M(x) + l 



-E 



Define /5(x) = e - 2 ^Q( x )/ x m We separate our analysis into two cases. 
Case 1: p(x) < p < 1. 

Note that for this range of values of p we have, by Lemma 3.3, that u(p) = — ^logpfl + o(l)), and by 

Lemma 5.1 (a.), that u{p{x)) ~ 2a < " 1 for all sufficiently large x. It follows that w(p,x) — 
min{u(p), lo^ 1 {x)/x\ = u(p). Also, by Lemma 6.1 we have that there exists a function tpi(x) \. as x — > oo 



such that 



\Z K (p, x) - A K (p,x)\ 

Oyfx{\ - p) 



E 



^ n=Af(a;)+l 



1(k > 2) 



1(« > 2) 



/ -i \ N ( x ) ^A„(u„) 

E ^^-e^»*(-7(p,*)) 

v r n=Af(a;)+l 
e f A pM»))$(- 7 (p, ^ > 2 ) + £ r p «(*.s) 1 ( CT Z < y/fll K (x)/y/x) 



(6.4) 
(6.5) 



:A p (u(p)) 



where 7(3;, p) = v / Mloga'/ ' + a^/xlogp/p 3 / 2 . Furthermore, by Lemma 6.2 we have that (6.4) and (6.5) are 
bounded by 

^(aOetM-to) ($(_ 7 ( P)X )) +^2(s)) +^ 2 (x)et A pW")) 

for some other y^O^) 4 0. Since by Lemma 5.2 we have that Z K (p, x) > Ce~> :KplyW ^ on p(x) < p < 1, it only 
remains to show that the term following (6.5) is o ^e^ Ap< - u ^ p ^ . First we notice that exact computation gives 

3 f Ap(-(p))$(_ 7 ( P; x )) i( K > 2 ) + E \p<^ z ~) 1 (aZ < ^lh k (x)/Vx~) 

:A p («(p)) 



\A p (u(p)) 



$(-7(p,x))l(« > 2) + e" 



f (logp+ -- ( '°f/' -A p (u(p)) 



3 ^A p («(p)) 



( 



y^i (g) CTy/alogp 



v p( u W)$( 7 (p,x)) 



log p+ - 



-A p (u(p)) 



K = 2, 
K > 2. 



Vph K (x) (Ty/X~\0gp 
<T\fx p?l 2 
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sup $ 

p(x)<p<l 



When k = 2 wc simply have 

y/^hHx) _ cry's log A = / ^//Za^Q:) 2gQ(x) \ 
cr^ M 3/2 ; ^ ^ + M V2 v ^ y | 

as x — > oo, since by Lemma 5.1 (a.) Q{x)/ui^[ (x) — > 0. When k > 2 note that 

~, , \ \ 1 6°sp+ g3 ty )2 -A P Kp))) -. 

sup $(7(/9,xj) e V 2tl / — 1 

<C sup $ -^-5 isn x |i ogp |3 

p(x)< P <i V CT / 



i)i 3 

(7 



c / VTHog^ 

< — = sup "P I 

V 1 0<t<2aQ(x)/^JIx \ 17 



as x — > oo. 

Case 2: < p < p(x) 



For this range of values of p we use Lemma 5.2 to obtain that Z K (p, x) > Cp(l — p) 1 e Q( x \ which together 
with Lemma 6.1 gives, for < p < p(x), 



\Z K (p,x)-A K (p,x)\ C(l-p)eQM \ {l- P ) ^ et A ^ 

— z^j) — - — - P — £ — ^> 2 ) 



e f 

n=M(x)+l 



+ E 



J^/xlogx u 

+e f ^p+ a - 1 ^ L ^ ( VE^M + ^iogp\ + eiAp(w(p , x)) 

\ a^fx p?l 2 J 

Let C = max{l,2/i 2 /er 2 } and p(x) = e"^ 1 ^)/^ and note that 

\Z K (p, x) — A K (p, x)\ _j Q ( s j _ s.A I) (a)r 1 (x)/a;) 

sup — -. r < sup Cp ' < yanogxe*' pv 1 w/ ' 

0<p<p(x) Z R (p,X) 0<p<p(x) I 

, e * lo gP+ ' 2 -^^ $ / CrV5|l0gpl / _ p?h K ( X ) 
\ p?l 2 \ £7 2 x|log|0| 



V^Uogpl 

< Ce Q(x) sup \ sc v / + -= re » logp ^ 

o< P <p( 2; ) [ \/a;|logp| J 
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which converges to zero as x — > oo since by Lemma 5.1, Q{x)/u 1 (x) — > and t/x/uj 1 1 (x) —> 0. For the 
range p(x) < p < p(x) we have, by Lemma 3.3, 

gup \Z K (p,x) - A K (p,x)\ 
p(x)<p<p(x) Z K (p,x) 

< sup CeW ly/Zlogxe* A 'M p)) 1(« > 2) + ^p+ ^X^' + e gM"00) j 

„ f-\ iE-log p+ ' r2x <'° . f p)2 (l+o(l)) 

< sup Ce yt J a;e" 

p(a:)<P<p(a:) 

= Cxe Q{x)+ » losp{x){1+o{1) ' ) = Cxe- Q{x){1+o{1) \ 
which also converges to zero as x — > oo since Q(x) > 21ogx by (2.2). This completes the proof. □ 



7. Numerical examples 



We conclude the paper with two examples comparing simulated values of P{Woo{p) > x) to the approxima- 
tions Z K (p,x) and A K (p,x) suggested by Theorems 3.2 and 3.4. For illustration purposes we also plot the 
heavy-tail and heavy-traffic approximations 

— — F(x) and expi--(l-p)!. 
l-p { p J 

The simulated values of P{W 00 {p) > x) were obtained using the conditional Monte Carlo algorithm from 
[2], and each point was estimated using 100,000 simulation runs. We point out that simulating heavy-tailed 
queues in heavy traffic is very difficult, and in particular, the simulated values of PiW^p) > x) for pairs 
(a;, p) in the region around the point where the queue's behavior transitions from the heavy traffic regime into 
the heavy tail regime, are highly unreliable. In terms of the approximations Z K (p, x) and A K (p, x) suggested 
in this paper, they tend to be sensitive to the mean and variance of the integrated tail distribution, p and 
er, respectively, so we suggest first scaling the queue in such a way that both parameters are small (of order 
one). We give two examples below, one in which the integrated tail distribution is lognormal and one where 
it is heavy-tailed Weibull; note that no M/G/l queue can have exactly Weibull integrated tail distribution, 
since its density is not monotone, but there are valid distributions (with decreasing densities) whose tail is 
asymptotically Weibull. For the lognormal(a, (3) example we used Q(x) = (log x — a) 2 /(2/3 2 ), which although 
an approximation to log F{x) works well in practice. 
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Fig 1. Lognormal(ct, 0) integrated tail with p = 0.9, a = 0, /3 = 1. 




Fig 2. Weibull{ct, j3) integrated tail with p = 0.9, a = 0.5, j3 = 0.22361. 
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C n (4.6) 
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